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Abstract. We explore the topology of configuration spaces of hard disks ex- 
perimentally, and show that several changes in the topology can already be 
observed with a small number of particles. The results illustrate a theorem 
of Baryshnikov, Bubenik, and Kahle that critical points correspond to config- 
urations of disks with balanced mechanical stresses, and suggest conjectures 
about the asymptotic topology as the number of disks tends to infinity. 

1. Introduction 

Configuration spaces of hard spheres are of intrinsic interest mathematically 
and physically, yet so far little seems known about their topology. As shown in 
this note, the dependence on the volume fraction (or equivalently on the sphere 
radius) of properties as fundamental as the connectivity of configuration space is 
already quite complicated even for a small number of spheres. More generally, we 
are concerned not only with the property of connectivity or with the number of 
connected components /3o, but also with the ith Betti number which denotes the 
dimension of ith homology, or colloquially "the number of i-dimcnsional holes pQ." 

One purpose of this article is to illustrate some of the complexity of the topol- 
ogy of configuration spaces of hard spheres, even in two dimensions (i.e. hard 
disks) when the number of disks is small. We illustrate a theorem of Baryshnikov, 
Bubenik, and Kahle that mechanically-balanced (or jammed) configurations act as 
critical points, indicating the only places where the topology can change 0. We 
always mean "critical point" in the mathematical sense of non-differentiable or sin- 
gular, rather than the statistical mechanical sense. Examples of such configurations 
may be seen in Figure [T] The combinatorial and geometric structure of the bond 
networks of these critical points is important in classifying their degeneracy and in 
computing their index. Furthermore, the bond networks suggest connections to the 
geometric structure approach in hard particle packing [3], as well as to tenscgrity 
and strut frameworks [4]. 

We show here that the topology of these configuration spaces can already be 
fairly complicated in two dimensions with a small number of particles. For example, 
with five disks in the unit square, the topology changes at least 20 times as the 
disk radius varies, and for radius 0.1686 < r < 0.1692 the configuration space has 
ft = 2176 "holes". 

The topology of configuration spaces has been studied in the context of motion 
planning for robots; see the surveys by Farber and Ghrist 0(3], and the textbook by 
LaValle [7J. One central idea in this field is that topological complexity, measured 
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Figure 1 . Index zero (bottom row) and index one (top row) crit- 
ical points, shown with their bond networks. 

in various ways, is a measure of the difficulty of constructing a motion planning 
algorithm. Many papers idealize robots as point particles in a suitable metric space 
— although it has been noted that the assumption that they have no volume is not 
physically realistic [8]. On the other hand, expanding the particles to have positive 
thickness greatly complicates the topology of the underlying configuration space — 
see for example the recent work of Deeley on thick particles on a metric graph [5] . 

As further motivation for this study, a set of hard disks inhabiting a bounded 
area may be regarded as a model thermodynamic system. From this standpoint, 
a topological property such as connectivity of the configuration space is a funda- 
mental concern, since this is required to satisfy the assumptions of equal a priori 
probabilities and ergodicity (e.g., Section VII of Torquato and Stillinger [3]) inher- 
ent to the statistical mechanics formalism. Despite the importance of this question, 
few results seem to be known for finite hard disk systems [TU] , and fewer still in the 
thermodynamic limit. It is known that if one can only move one disk at a time, the 
configuration space for large numbers of disks can be disconnected even for density 
arbitrarily close to zero [TTJ [T2] . Even if collective motions of disks are allowed, 
the configuration space may be disconnected for density bounded away from maxi- 
mal — for example, the reinforced Kagome lattice is strictly jammed and only has 
density Ri 0.68 compared to the close packing density of ~ 0.91 [T5] , 

Perhaps one of the most notable properties of the hard disk system in the con- 
text of thermodynamics is the existence of a solid-to-liquid phase transition [14] 
with the variation of the packing fraction, or equivalently of the disk radius. This 
is particularly suggestive in the context of a number of recent papers that ex- 
plore the hypothesis that statistical-mechanical phase transitions are intimately 
related to changes in the topology of surfaces of constant potential (i.e. equipo- 
tential submanifolds) in the configuration space [15J [16j H3 HB]- Specifically, a 
theorem of Franzosi, Pettini, and Spinelli jTHK^O] states that for rather benign con- 
ditions (smooth, finite-range, confining potentials), the low-order derivatives of the 
Helmholtz free energy cannot display a discontinuity unless there is simultaneously 
a change in the topology of the equipotential submanifolds with the variation of a 
relevant parameter, such as the disk radius. While this indicates that a topologi- 
cal change is necessary, a conjecture known as the Topological Hypothesis roughly 
states that a major topological change (a proposal for specific conditions is given by 
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Kastner, Schreiber and Schnetz [2"H 122] ) is sufficient to drive the phase transition 
as well. This immediately raises the question of the topological signature of the 
solid-to-liquid phase transition for the hard disk system. 

Historically, the accumulation of certain indirect evidence encouraged the devel- 
opment of the Topological Hypothesis. For example, the magnitude of fluctuations 
in the curvature of equipotential submanifolds was observed to show singular be- 
havior near a phase transition |23j . This singular behavior seemed to be insensitive 
to the choice of metric [21] , pointing to a more fundamental change in the topology 
of the accessible portion of configuration space at the phase transition. Further 
support was provided by the existence of major changes in the topology of the 
equipotential submanifolds accompanying phase transitions in certain exactly solv- 
able systems [25J 123 HE1 HZ] ■ Strikingly, the presence of these changes is known to 
persist for small systems well away from the thermodynamic limit 25, 2*f\. That is, 
a major topology change in the equipotential submanifolds of a system with just 
a few hard disks may indicate the presence of a topology change that drives the 
phase transition of the corresponding system in the thermodynamic limit. 

There is a persistent difficulty faced by investigations of the Topological Hypoth- 
esis in the physics literature, namely, that the exactly solvable systems mentioned 
above do not satisfy the conditions for the theorem of Franzosi, Pettini and Spinclli. 
Specifically, they violate the requirements of short-range interactions [S5J [551 HZj 
and confining potentials [16] . This stands in sharp contrast to the short-range, 
hard-core repulsive potential characteristic of hard disks systems, meaning that a 
complete characterization of the configuration space topology of our system will 
begin to address a need that is openly acknowledged in the literature [T7] . 

We begin this study by reconstructing the dendrogram, i.e. a graph which en- 
codes the appearance of path components and the connections between them as 
the disk radius (or volume fraction) varies. More generally, for a suitably defined 
energy function, the dendrogram is defined to be the graph with nodes given by the 
basins of attraction in the energy landscape and edges given by the saddle points 
that connect the basins as the energy increases. These correspond in the physics lit- 
erature to collectively-jammed configurations and to the transitions between them 
(see Torquato & Stillinger [3]), subjects of continuing interest to the physics com- 
munity. The name "dendrogram" comes from hierarchical clustering in statistics. 
We find the saddle points using the "nudged elastic band" method [25] ■ From the 
point of view of Morse theory, the index of the basins of attraction is zero since 
there is no collective motion that decreases the energy, and the index of the saddle 
points is one since there is one dimension of collective motion that decreases the 
energy. Examples of these varieties of critical points appear in Figure [l] 

As alluded to above, these features of the configuration space correspond in the 
Morse-theory literature to critical points of index zero and index one, respectively 
[29] . Along with the notion of a critical index comes the notion of higher index 
critical points, allowing for the construction of a more complete picture of the 
configuration space and the related energy landscape. This is precisely the field of 
Morse theory, which provides a framework to pass from the populations of critical 
points with specified indices to a reconstruction of the configuration space, and 
allows us to answer many more questions than simply the energy at which the saddle 
point connecting jammed configurations is accessible. For instance, we deduce 
certain facts about the Betti numbers /?/. that indicate the number of fc-dimensional 
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holes or the dimension of kth homology, and observe suggestive patterns about the 
asymptotic topology of configuration spaces for large numbers of disks. 

There is good reason to introduce Morse theory as a fundamental instrument 
for the study of configuration spaces. Various flavors of Morse theory sit at the 
heart of Franzosi, Pettini, and Spinelli's proof [19] in the context of the Topological 
Hypothesis, Yao et al.'s exploration of low-density states in biomolecular folding 
pathways [30], and Deeley's study of configuration spaces of thick particles on a 
metric graph [9]. Section [3] will briefly overview Morse theory, though most of the 
technical mathematical details, and particularly the proof of the main theorem we 
quote, are left to another article 2 . Our primary focus in this article is on our 
numerical experiments to identify critical points of various indices and connecting 
the Morse theoretic viewpoint of configurations spaces with related concepts in the 
physics literature. 

2. Configuration spaces 

Let C(n, r) denote the configuration space of n nonoverlapping disks of radius 
r in the unit square [0, l] 2 . If we indicate the coordinates of the center of the ith 
disk by X{, we can write this as 

C(n,r) = {(xi,X2,...,x n ) | Xi € [r, 1 - r] 2 , d(x t , xj) > 2r}. 

Since C(n,r) C R 2n , it naturally inherits a subspace topology and a metric from 
the Euclidean space of the appropriate dimension. The configuration space of hard 
disks in the limit of large r is the empty space, since the disks do not fit in the 
unit square. Meanwhile, the configuration space of n labelled points in the plane 
(or more generally on a manifold) is well studied |31L I32j , and is equivalent to the 
configuration space of hard disks in the limit of small r. This leaves the study of 
the topology of C(n, r) for intermediate r, which is much more complicated than 
for the two cases above. 

Notice that if r' < r, then there exists a natural inclusion map 

C(n,r) ^ C(n,r') 

given by shrinking the disks with fixed centers. In other words, as the disk radius 
shrinks, the space of allowable configurations grows. 

As mentioned in Section [Tj a basic question is whether C(n, r) is connected for 
a given value of r, and if not, how many components there are. A more subtle 
description involves the construction of the dendrogram, or the graph that encodes 
the appearance of persistent components and the connections between them as the 
disk radius varies. Equivalently, the dendrogram is a means to visualize distinct 
basins of attraction in the energy landscape and the saddle points between them, 
not the properties of the configuration space at higher energies. That would require 
identifying the critical points of all indices and using Morse theory to reconstruct 
the remaining features of the configuration space. 

We give a brief discussion of Morse theory to establish a unifying set of language 
and ideas. Following that, we apply each of these approaches to the specific case 
of the configuration space of five disks in the unit square. 

3. Morse theory 

The basic idea of Morse theory is to filter a topological space X by the sub- 
level sets oo, r]) of a suitable function / : X — > K [551 IS3]- Suppose that 
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the function / assigns an energy E to every point in the configuration space, and 
that this energy increases as the radius of the disks decreases. The sub-level set 
associated with a particular energy Eg is just the set of accessible states, or the 
set of points in the configuration space where the energy is less than or equal to 
Eq . The existence of the inclusion map in the preceding section is equivalent to the 
statement that as Eq increases, the set of accessible states increases. Provided that 
the function / is well-behaved, the topology of the sub-level sets will change only a 
small number of times as r decreases. One of the tenets of traditional Morse theory 
is that these changes occur precisely at the critical points where V/(x) = 0, or at 
the mechanical equilibria of the system. Using techniques similar to those in "Min- 
type" Morse theory [34] and stratified Morse theory 35j, Baryshnikov, Bubenik, 
and Kahle recently proved the following theorem [2J. 

Theorem 3.1. If there are no "balanced" configurations of disks of radius r with 
r € [a,b] then C(n,a) is homotopy equivalent to C(n, b). 

The definition of a "mechanically-balanced" configuration is as follows. For a 
configuration of disks, define the bond network to have vertices at disk centers, and 
edges e.ij connecting pairs of tangent disks and also disks to the boundary 

of the region. We consider the edge ey to be a vector, pointing from i to j, so in 
particular = — e^. 

Definition 3.2. A configuration is mechanically-balanced if non-negative weights 
Wij can be assigned to the edges of the bond network so that at least one Wij is 
positive, and for every fixed vertex i, 



where the sum is over all vertices j adjacent to i. 

The point here is that C(n,r) is the sub-level sets of some function — but not 
a smooth function. Smooth Morse theory and Morse-Bott theory don't quite work 
in this case, but something similar to min-type Morse theory |34| does work. We 
note that this is similar in spirit to the PL Morse-Bott Theory in Deeley's recent 
work (9]. 

For a smooth Morse function the index of a non-degenerate critical point is the 
number of negative eigenvalues of the Hessian. For a min-type Morse function 
(such as for hard disks), you can compute the index of critical points with linear 
programming, i.e. the index is the dimension of a certain cone. In either form of 
Morse theory, the fundamental theorem states that every time you have an index k 
critical point, the homotopy type changes by gluing in another /c-dimensional cell, 
and importantly — that there are not any other changes in the topology except 
when there are critical points. 



It is convenient for our methods to first slightly "soften" the disks. We recover 
their hardness shortly. In general, on the configuration space of n distinct points 
in a region R we can define a smooth energy function E by 




4. Part I: Dendrogram 
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where h ^> is a fixed "hardness" parameter and dR denotes the boundary of R. 
For our experiments we set h = 50, n = 5, and R = [0, l] 2 . 

Taking uniform random points in [0, l] 10 and flowing in the direction of steepest 
descent along — V-B, we find the three qualitatively different types of local minima 
illustrated in the bottom row of Figure [TJ We believe this to be a complete list of 
all "types" of local minima (see discussion below) for five disks in a box. 

Notice that a point in [0, l] 10 gives the coordinates of five distinguishable disks. 
Since exchanging the labels of the disks does not change the energy of the configura- 
tion, this permutation symmetry entails that every critical point be symmetrically 
equivalent to 120 points in the configuration space (including itself). (To mathe- 
maticians, the labeled configuration space admits a free action of the symmetric 
group, and the quotient is the unlabeled configuration space.) This number of in- 
stances of each type is further increased by the additional configurations generated 
by applying the symmetries of the unit square. For instance, the class of configura- 
tions symmetrically equivalent to the critical point at a radius of 0.2071 in Figure 
[2] contains 120 points, while the equivalence class containing the critical point at a 
radius of 0.1964 contains 480 points, and the equivalence class containing the criti- 
cal point at a radius of 0.1693 contains 960 points. So by a "type" of critical point, 
we mean an equivalence class of critical points under both types of symmetries. 

Once a local minimum of E is found, the preliminary radius r of the configuration 
is the maximum possible radius such that that no disks overlap or extend outside 
the unit square. The "bond network" of the configuration is found by connecting 
centers of disks that are within (2 + e)r of each other, or within (1 + e)r of the 
boundary. Here e > is a fixed small tolerance parameter, usually around 0.001. 
With the bond network established, we effectively harden the disks by searching 
for an exact solution to the set of equations given by varying r and setting e = 0. 
In some cases this gives a unique solution, and establishes the final values of the 
radius for the three red configurations along the bottom of Figure [2j 

We then consider pairs of local minima which are close in the configuration space 
for the induced metric and find low energy paths between them using the "nudged 
elastic band" method [351 [2S] • The highest energy points along a low energy path 
should be close to a saddle point. Indeed, by establishing the bond networks for 
all of the highest energy points along low energy paths and applying the hardening 
procedure described above, we identify the saddle points depicted in orange along 
the bottom of Figure [2j This allows us to construct the dendrogram for our system, 
though a full graphical representation of this graph would be quite complicated. We 
present simplified versions of some of the salient features in Figures [2j [4] and [5] 

A precise value of the radius is important to be able to order the critical points 
by the radius of the disks when they appear. Evidence that a rough numerical 
approach is not sufficient is provided by the two critical points of index one in 
Figure [3j where the radius differs only in the fourth decimal place. This is one 
reason why the more elaborate approach described above is necessary. 

Figure [2] shows the variation in the number of connected components as a func- 
tion of the radius. The total number of connected components for a given radius is 
shown at the top of the figure. The central part of the figure collects the connected 
components into classes of symmetrically equivalent elements, and labels the classes 
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Figure 2. (Color online.) A bird's-eye view of the dendrogram for 
labeled disks. New components appear at radius 0.2071, 0.1964, 
and 0.1705, and components merge together at radius 0.1942, 
0.1871, and 0.1693. 




0.1693 0.1692 

FIGURE 3. Two saddle points with nearly the same radius 



by the number of elements they contain. The bottom of the figure shows the con- 
figurations of the disks at the relevant critical points, and provides the values of 
the radius when these configurations appear. 

For r > 0.2071, the disks do not fit inside the unit square and the configuration 
space is empty. The first allowed configurations occur at a radius of 0.2071, when 
120 connected components appear simultaneously. All of these belong to the same 
class, and correspond to different labelings of the disks in the first red configuration. 
As the radius decreases, a class of 480 more connected components appear at radius 
0.1964, corresponding to different labelings and arrangements of the disks in the 
second red configuration. The 480 components of this class connect in groups of 
20 with the emergence of the first saddle point at radius 0.1942, leaving a class of 
24 embedded circles. One of these circles is shown in Figure [4| Every arc segment 
connecting the jammed configurations in Figure [4] is a low energy path containing a 
saddle point of the type appearing at radius 0.1942. The appearance of the second 
saddle point at radius 0.1871 connects all of the remaining components. The space 
remains connected for all smaller values of the radius, with the exception of the 
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FIGURE 4. (Color online.) For radius 0.1871 < r < 0.1942, 
each connected component of the configuration space is homotopy 
equivalent to a circle. One such embedded circle and 20 index-0 
critical points is shown here. Not shown are the 20 transitionary 
index- 1 critical points between them. 



^13 ET? 




Figure 5. (Color online.) The dendrogram for unlabeled disks. 
The numbers indicate the order in which critical points appear as 
the disks shrink — new components appear at 1, 2, and 5, and 
components merge at 3, 4, and 6. The radius where each vertex 
and edge appear is given in Figure [HJ 



appearance of 480 more components at a radius of 0.1705, though these connect up 
with the main component soon afterward. 
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An alternate representation of the dendrogram is given in Figure [5] Since the 
purpose of this figure is to allow the visualization of the paths connecting energy 
minima, we consider only unlabeled disks for the sake of simplicity. The labels of the 
critical points indicate the order of appearance as the radius of the disks decreases, 
where a set of critical points with the same labels constitutes an equivalence class. 
Similar to the description in Figure [2] the first allowed configuration to appear 
is the red critical point in the center of the figure. As the radius decreases this 
single configuration is joined by the four critical points arrayed around the central 
circle, though they all remain isolated until the first saddle points appear. At this 
interval the configuration space is given by the isolated configuration in the center 
and a circle that is roughly analogous to that in Figure |4j but with unlabeled disks. 
The central configuration and the circle become connected with the appearance of 
the second saddle point, and the space remains a single component until the final 
energy minima appear at the outer extremities of the figure. The third saddle point 
once more connects the space into a single component. 

Figures [2] and [5] do not show all of the critical points of index one for five disks 
in the unit square. Indeed, we identified five types of critical points of index one, 
though only three of these change the number of connected components in the 
configuration space. Figures [2] and [5] show precisely these three types of critical 
points, since the remaining two do not directly affect the connectivity of the space. 

5. Part II: Higher index critical points 

While the above section is concerned with the critical points of index zero and 
index one and the relation of these to the dendrogram, the current section is con- 
cerned with more general features of the configuration space. As mentioned in 
Section [T] this requires the identification of the critical points of higher index. For 
this purpose, we work with the same smooth function E defined above. Using 
the MuPAD computer algebra system included in the Symbolic Math Toolbox for 
MATLAB, we symbolically differentiate E and set 

F = \\VEf. 

Since VE = at the critical points of E, F = at the critical points of E as well 
and is strictly greater than zero everywhere else. Hence, we may find critical points 
of E of any index by flowing down the function F to one of the zeros. 

More specifically, our search begins from a random point in [0, l] 10 . We symbol- 
ically calculate the gradient of F with MuPAD and apply the conjugate gradient 
method with the Polak-Ribiere formula, using the secant method to find the min- 
imum of F in the search direction on every iteration. This gives much more rapid 
convergence to the local minima of F than the steepest descent algorithm. 

Once the conjugate gradient algorithm converges, the preliminary radius r of 
the configuration is the maximum possible radius such that that no disks overlap 
or extend outside the unit square. The bond network is constructed just as for the 
calculation of the dendrogram, though the tolerance parameter e is on the order of 
0.01. We recover the configuration of the critical point in the hard disk system by 
modeling the bonds in the bond network as springs. Let s be the average of the 
bond lengths from the disks to the boundary and half the bond lengths connecting 
pairs of disks. Set the equilibrium length of the springs from the disks to the 
boundary to be s, and the equilibrium length of the springs connecting pairs of 
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Figure 6. From left to right, critical submanifolds of dimensions 
1, 2, 3, and 5. 



disks to be 2s. We minimize the potential energy of the spring network using the 
conjugate gradient method and update s on every iteration. This is repeated while 
reducing the spring constant from a small positive number to zero, allowing us to 
find the radius of the disks at a critical point of the hard disk system to more than 
five significant digits. 

The index of a critical point is found by symbolically calculating the Hessian of 
E, evaluating this expression for the given disk locations, and numerically finding 
the number of negative eigenvalues. Provided that the Hessian is well-conditioned, 
the index of a critical point calculated for the smooth function E should be the 
same as for the hard disk system. Some discretion should be used when finding the 
critical index this way though, since the Hessian is not always well-conditioned. In 
these cases, the index is found as the dimension of a certain cone calculated from 
the bond network. 

We found only 30 types of critical points and critical submanifolds after repeating 
this process on the order of 10 5 times. By "critical submanifold" we mean a critical 
point with several degrees of freedom, where the number of degrees of freedom is the 
dimension of the submanifold. (In the strict mathematical sense, these might not 
be manifolds or even manifolds with boundary.) Specifically, the presence of each 
"rattler" contributes two dimensions to the critical submanifold, while a row of disks 
across the square adds one dimension. Several examples of these appear in Figure 
|6j along with the corresponding dimensions. Notice that the 5-dimensional sub- 
manifold on the right is not contractible, as the others are — this one is homotopy 
equivalent to a circle. 

Of these 30 types of critical points and submanifolds, we consider the 10 in Figure 
[7] to be degenerate. Our criterion is if weights Wij may be assigned to the edges of 
the bond network of a balanced configuration so that at least one Wij is positive 
and at least one tUy vanishes, then the critical point or submanifold is degenerate. 
Computing the local topology near several critical points that satisfy this condition 
has convinced us that the topology does not change at these critical points, but we 
do not yet state this as a formal theorem. We intend to explore the question of 
degeneracy in more mathematical detail in a future paper. 

The remaining 20 non-degenerate critical points and submanifolds appear in 
Figure |8j ordered horizontally by the radius of the disks and vertically by index. 

This figure reveals that higher index critical points generally coincide with smaller 
radius. For instance, all of the index-0 and index-1 critical points occur at larger 
radius than any of the higher-index critical points. The fundamental theorem in 
Morse theory states that each critical point of i contributes an i-dimensional cell 
via an attaching map, up to homotopy equivalence (smooth deformation) |35[ 129] . 
This indicates that for radius r > 0.1692, the whole configuration space is homotopy 
equivalent to a 1-dimcnsional cell complex, i.e. a graph. 
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Figure 7. Degenerate types 
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Figure 8. (Color online.) Non-degenerate types, in rows by index. 
Index critical points are in the bottom row, all the way up to 
index 4 in the top row. 



The index-0 critical points correspond to vertices of this graph, and the index- 1 
critical points correspond to edges. Since we have already computed the number of 
connected components Pq in the dendrogram above, the Euler relation 

v - e = p ~ Pi ■ 

allows us to compute Pi — the results in Figure [9] While the topology changes 
from the two critical points with the same radius r w 0.1705 are simultaneous, 
their indices are distinct and their contribution to the Bctti numbers are computed 
separately. The calculation does not continue for r < 0.1686 since we only have 
bounds on the Betti numbers after this point — with the attachment of 2-cells, the 
Euler relation 

v - e + f = & - Pi + A> 
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Figure 9. (Color online.) Betti numbers for the first several 
changes in topology. 



alone is not sufficient to find fi\. 

From Figure [9] we see that for 0.1686 < r < 0.1692 the configuration space is 
homotopy equivalent to a bouquet of 2761 circles. Since we know that no more 
1-cells are attached for r < 0.1692, we have that j3\ is monotone decreasing as the 
radius decreases from this point on. Meanwhile, for small enough radius, j3\ = 10 
since this is the Betti number for the configuration space of 5 points in the plane. 
Therefore, for five disks in the unit square we know that fix is a unimodal function 
of the radius and that its peak value is j3i — 2761. 

The non-degenerate critical points are shown sorted by radius (and in rows by 



index) in Figure 10 In this histogram, it is seen that all the critical points have 
radius between 0.100 and 0.220 — however half of them occur between radius 0.160 
and 0.175. This suggests the intriguing possibility that a large number of critical 
points are concentrated in a small interval as the number of particles n — > oo, 
potentially the type of situation required by the Topological Hypothesis to drive a 
phase transition. 

We do not know with certainty that we have a complete list of all types of critical 
points for 5 disks in a box — for example, it is easily possible that some points 
have very small basins of attraction for the algorithms we used, and that we also 
failed to find them by hand. If so, it would change the calculations we have made 
— nevertheless, we expect that we have most of them and one or two missing types 
would likely not significantly alter the qualitative behavior we observe. 



6. Comments 

It is probably within computational reach to implement experiments like these 
for greater numbers of disks, spheres in higher dimensions, and in other types of 
bounded regions. Flat tori and sphere are especially attractive boundaryless 2- 
dimensional cases. Even if it is difficult to find a complete list of critical points 
for large numbers of disks, it would still be interesting to sample a large selection 
of them. It might be possible, for example, to find a geometrically feasible bond 
network and then increase the disk radius, preserving tangency at all time, until 
a critical point is found. Compiling statistics on radius and index might already 
give interesting information, as in Figure |10| New techniques for exhausting index- 
jammed configurations, and over all possible lattices, are described in [37], and 
these methods could likely be adapted to higher index points. 

A promising area for future study is the asymptotic topology of C(n, r) as n — > 
oo. It is generally not even known how large one can make r = r(n) and maintain 
the property that C(n, r) is connected as n — >■ oo. For ergodicity in a convex region, 
it is straightforward but important to ensure that the condition r < C/n for some 
region-dependent constant C implies that C(n,r) is connected [35] [3J5]. However, 
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Figure 10. Histogram of the nineteen types of non-degenerate 
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spheres satisfying this condition are quite small, taking up only length rather than 
volume in the limit, and it may be that some configuration spaces are connected 
for much larger radius. 

The question of threshold radius for connectivity depends on the bounding region 
- the case of a circle implies that the C/n bound is the best possible when the 
region is not constrained. It has been conjectured that for a flat torus though, the 
least dense collectively jammed configuration is the (reinforced) Kagome lattice [T5] , 
If an analogue of this conjecture holds for index zero points and index one critical 
points, for example, then r = 0(n~ 1 / d ) would ensure connectivity of C(n,r). 

A more general question is the asymptotic behavior of the Betti numbers when 
the radius is set to r = r(n) and n — > oo. We can guess but not yet prove several 
of the qualitative features of the asymptotic topology of C(n, r) for this situation. 
For example, we expect /3q to be roughly decreasing with radius, and each Betti 
number j3k with k > 1 to be roughly unimodal in the disk radius. In support of this 
conjecture, we display the Betti numbers we computed for four disks in a square 
in Figure [TT] — in this case there are fewer changes in the topology, the index of a 
critical point is monotone decreasing in the radius, and we were able to compute 
all the Betti numbers. 
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